The ac magnetic response of mesoscopic type II superconductors 



Alexander D. Hernandez 1,2 and Daniel Dommguez 2 

1 Laboratorio de Superconductividad, Facultad de Fisica-IMRE, 
Universidad de la Habana, 10400, Ciudad Habana, Cuba. 
2 Centro Atomico Bariloche and Instituto Balseiro, 
8400 San Carlos de Bariloche, Rio Negro, Argentina 

The response of mesoscopic superconductors to an ac magnetic field is numerically investigated on 
the basis of the time-dependent Ginzburg-Landau equations (TDGL). We study the dependence with 
frequency u> and dc magnetic field Hdc of the linear ac susceptibility x(Hdc,u) in square samples with 
dimensions of the order of the London penetration depth. At Hdc = the behavior of x as a function 
of uj agrees very well with the two fluid model, and the imaginary part of the ac susceptibility, x" (w), 
shows a dissipative a maximum at the frequency v — c 2 /(4na A 2 ). In the presence of a magnetic 
field a second dissipation maximum appears at a frequency cu p <C vq. The most interesting behavior 
of mesoscopic superconductors can be observed in the x(Hdc) curves obtained at a fixed frequency. 
At a fixed number of vortices, x" (Hdc) continuously increases with increasing Hdc- We observe that 
the dissipation reaches a maximum for magnetic fields right below the vortex penetration fields. 
Then, after each vortex penetration event, there is a sudden suppression of the ac losses, showing 
discontinuities in x" (Hdc) at several values of Hdc- We show that these discontinuities are typical 
of the mesoscopic scale and disappear in macroscopic samples, which have a continuos behavior of 
x(Hdc)- We argue that these discontinuities in x(Hdc) are due to the effect of nascent vortices which 
cause a large variation of the amplitude of the order parameter near the surface before the entrance 
of vortices. 
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I. INTRODUCTION 



The response of superconductors to a: 
field has been of interest fehi^^R^tiflpS 
ticularly in the last years.ffiffl'fflElElB 
crowave surface impedance Z s = R s — iX s and the ac 
magnetic suscepUhil^^^,-^^^'' have been ex- 
tensively studied.HffliyBQ&BByiae-y Interesting be- 
havior has been found for different values of the fre- 
quency, the magnetic Ji^ld, the pinning force and the 
thermal fluctuations jStjaLJ whereas linear and nonlinear 
response appear depending on the strength of the ac 
signal.Q The main interest so far has been on the elec- 
trodynamics of macroscopic samples containing a large 
number of vortices. In this case, it is possible to use phe- 
nomenological models to describe thamacroscopic behav- 
ior of the superconducting samples .oBQiII 

Recently there has been an increasing interest 

scale@l3SS0H3SSiS@BH a regim^in 
which a small number of vortices are confined in a 
small sample. The behavior of mesoscopic supercon- 
ductors is different from the behavior of bulk saixt. 
pies. In mesoscopic superconductors surface effectsEHl 
are very important since the interaction of vortices 
with the surface currents is large.E3 The magnetic 
properties_j3irp-Dgly depend on the sample sizes and 
geometryOOEf^i^pmi^^to are the the 

^idlidbl^yy^yi^^ m e S oscopic 

In particular 



mes pscc vpic dist 

slablijO and the mesoscopic square.! 
the mesoseepic samples can develop Abrikosov multivor- 
tex statescil and depending on the size of the sample it is 



possible to observe first or second order transitions Exlha 
electric charge of vortices£3 a paramagnetic effect]i3'E£l 
and the surface barriercj£3 on a mesoscopic scale have 
also been studied. One interesting characteristic of the 
magnetic properties of mesoscopic superconductors is the 
behavior of the dc magnetization curves. In a mesoscopic 
scale, the vortices that are inside the sample induce a re- 
inforcement of the surface barrier at fields greater than 
the first penetration field H p !Ei These new barriers allow 
for the existence of metastable states of constant vorticity 
as a function of magnetic field. Each metastable state be- 
comes unstable at the z-th penetration field H p i in which 
vortices enter the sample and the magnetization has a 
discontinuous Jimrp. These jumps have been observed 
numepcallyo'E2fE3 and experimentally in mesoscopic Al 
disksEJ 

In the last years numerical simulations of the time de- 
pendent Ginzburg Landau (TDGL) equations have been 
an important tool in the study rfrthe stotifl ajuLxfoc 
properties of superconductors. ESEiJoli 
The TDGjL.xn.odel has been used to study 
dynamicsEH, the magnetic responseEj'l 
characteristics of sytaerconducting samples £3 In particu- 
lar, Enomoto et alBB using the TDGL equations studied 
the temperature dependence of the ac susceptibility at a 
fixed frequency in the absence of a dc magnetic field for 
a large sample. For increasing temperature they found 
that x' exhibits a step-like change from a negative value 
(x' = — 1/(4tt) to zero, while x" initially rises from zero, 
goes through a maximum, and then returns to a small 
value near T c , which is qualitatively consistent with the 
overall behavior observed in experiments. 

In this paper we perform numerical simulations of 
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the time dependent Ginzburg-Landau equation to study 
the ac magnetic response of type II superconductors on 
a mesoscopic scale. We study the frequency and Hd c 
dependence of the linear ac susceptibility x(Hd c ,u>) — 
x'(HdctU) + *x" (Hdc,u) of square samples with dimen- 
sions of the order of the London penetration depth Al. 
Our resuls were obtained in the absence of bulk pinning 
and in the linear regime. 



The paper is organized as follows. In Sec. II we review 
the known results for the ac magnetic response of macro- 
scopic superconductors, in order to later compare them 
with our results in mesoscopic samples. In particular, 
the two fluid model for the Meissner state and the Cof- 
fey and Clem model for the mixed state are reviewed. 
In Sec. Ill we describe the TDGL equations used in the 
simulations. In Sec. IV we present our numerical results 
for the ac magnetic response of mesoscopic square sam- 
ples. Sec. IV. A is devoted to the study of the frequency 
dependence of x' an d x" ■ We find that for Hd c = the 
behavior of the sample is well described by the two fluid 
model. At u> = v Q = c 2 / (AiraX 2 ) a maximum in the imag- 
inary part of x(co, 0) appears. In the Meissner state for 
Hdc ~ Hp, another maximum appears in x"(w). These 
behavior, characterized by the presence of two dissipa- 
tion maxima is also observed in the mixed state. At the 
same time, for increasing frequency x'( w j Hd c ) changes 
from London screening to a perfect diamagnetic state at 
high frequency. This transition is due to a decrease of 
the ac penetration length. In Sec. IV. B we study the Hd c 
dependence of x' an< 4 x" a t a fixed frequency. Depend- 
ing on the sample size, we find well defined mesoscopic 
and macroscopic behaviors. For mesoscopic samples we 
observe that, in the magnetic field ranges in which the 
number of vortices is constant, x" (Hdc) continuously in- 
creases for increasing Hd c - At the z-th vortex penetration 
field H p i, the entrance of vortices produces a consider- 
able suppression of the ac losses and x" (Hdc) decreases 
with a discontinuous jump. On the contrary, samples of 
macroscopic sizes show a continuous behavior in x(Hdc), 
with x" (Hdc) monotonically increasing with Hd c - In Sec. 
V we study the time evolution of the order parameter 
in mesoscopic superconductors. We show that at high 
frequencies and in mesoscopic samples, the vortices are 
fixed and only play a secondary role in the ac response 
of the sample. We show that the main dissipation mech- 
anism in mesoscopic superconductors is due to the effect 
of "nascent vortices" . They cause large variations of the 
amplitude of the order parameter at the boundary of the 
sample before the entrance of vortices. Finally in Sec. VI 
we give a summary of our results and conclusions. 



II. THE AC RESPONSE OF MACROSCOPIC 
SAMPLES: REVIEW OF KNOWN RESULTS 

A. Two fluid model and Coffey and Clem model 

When an ac magnetic field is applied H (t) = Hd c + 
h ac cos(ujt), an effective complex penetration depth X(u>) 
can be defined assuming that the fluctuating field inside 
a scmiinfinitc sample for x > has the form 



SH = h nr e- x / x ^e' luJt 



(1) 



or from a generalization of London's expression for the 
current, 



V x 3(lo) 



(2) 



47rA 2 (w) 

It can be related to the frequency dependent conductiv 
ity as A(w) = J j^J) 



In the case of a normal metal 



with conductivity <r„, the complex penetration depth is 



directly related to the skin depth S n = 



as \ n — 



(1 + i)5„/2. The complex permeability p-(w), or the com- 
plex susceptibility x(uj) — x'(w)+«x"( w ) = [m( w )~1]/47t, 
depend both on A(w) and also on the specific shape of the 
sample, since the ac penetration of magnetic fields is a 
surface phenomenon. In the case of a square sample of 
size LxLan approximate expression for the permeability 
in the limit A <C L is 



4A 
T 



1 - 



A 



(3) 



Exact expressions for jl for squares, cylinders and slabs 
can be found for example in Ref.^. Sometimes the com- 
plex surface impedance Z s is used, it can be defined when 
A < L as Z s = R s - iX s = iu^X(uj) . 

In the Meissner state the ac response -has been usually 
described with the two fluid model Jf-reHl In the mixed 
state, the magnetic flux enters in the form of quantized 
vortices. Therefore the dynamics of the vortices has to 
be included in a description of the magnetic ac response 
of the mixed state.QBtrEl In particular, Coffey and Clem 
have extended the two fluid model by including the. equa- 
tion of motion of vortices for small [displacements. □ Simi- 
lar results were obtained by Brandtu and by van der Beek 
et alB All these models assume that there is a large num- 
ber of vortices in the sample and a continuous description 
of the vortex lattice is used (valid for B ^> H c \). 

The two fluid model consists in writing the total cur- 
rent as the sum of the supercurrent and the normal cur- 
rent, 



with 



(T„E 



cr n dA 

' c at 



(4) 



(5) 
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and the supercurrent given by the London model 



C. Vortices without pinning 



V x J s 



(B 



(6) 



Here Xl is the static London penetration depth and b„ 
is the local vortex magnetic field. In the absence of vor- 
tices, h v = 0, the two fluid model has the characteristic 
time to — 47rA 2 <7„/c 2 for the transformation of supercur- 
rents into normal currents, since J n = todJ s /dt. The 
Coffey and Clem model includes the equation of motion 
of vortices for small displacements u(x, t) from their equi- 
librium positions, 



r/yii + K p u = - J x $ q: 



(7) 



where rj v is the viscous drag coefficient, t\ip IS the restoring 
force constant (Labusch parameter) of a pinning poten- 
tial well and 6t the local vortex direction. Here we have 
neglected thermal fluctuation effects. The local vortex 
magnetic field depgnds on the vortex displacements as 
b„ = V x (u x B).Q 

In the case of fields H || h ac and parallel to the surface, 
using Jilqs.^),^-^ the effective complex penetration 
depthta can be obtained as 



A 2 = A|- 



1 — iujto 



. 2X 2 L 

1 S 2 



(8) 

with the Campbell penetrationd depth X c = B^q/4ttk p , 
the flux- flow time scale tff = 4-7rA 2 <j///c 2 , the flux- 
flow conductivity <Jff = c 2 7] v /B$>o and the flux-flow skin 
depth $ff = c ' 2 /^ 7rLL " T ff- Since er// ~ a n H C 2/B, we have 
Off > a n and tff > to. Therefore the highest charac- 
teristic frequency is v$ — 1/to, which for conventional 
superconductors is about vq ~ 10 — 100 GHz.c3 



B. Meissner state 



In the absence of vortices the complex penetration 
depth is simply A 2 = A|/(l — iujt ). This leads to a 
dissipation peak in x" when ujto re 1. For large frequen- 
cies, cut 3> 1 (i.e. S n <C Xl), the system behaves as a 
normal metal (normal state skin depth effect) with: 



A' = A" re X, 



1 



2uito 



(9) 



Therefore dissipation goes as x" ^ X" / L ~ lu~ x / 2 for 
large ui. For low frequencies, cuto <C 1 (i.e. S n S> Xl), the 
system is dominated by the Meissner effect and 

1 



X'kXl, A"«Ai^ (10) 
and therefore dissipation goes as x" K X" /L ~ uj for low 



In the presence of vortices, the interesting frequency 
range is coto <C 1 since for frequencies Lut ^> 1 the system 
always behaves as in Eq.(^J), corresponding to the normal 
state skin depth effect. 

Let us first discuss the case when there is no bulk- 
pinning, k p = (Ac = oo). For frequencies such that 
utff ^ 1, the system behaves similarly to the Meiss- 
ner state, with a dissipation peak at ujto re 1, as de- 
scribed previously. For frequencies such that u>tff <C 1 
the system is dominated by the "flux flow skin depth ef- 
fect" due to the flux flow conductivity 07/ . In this limit 
A 2 re i8 2 ff/2 — iX 2 L /ujtff, and therefore the susceptibility 

for low frequencies should diverge as x ~ w -1 ' 2 . This 
divergence is cut-off by the finite system size. At low fre- 
quencies the real part of the effective penetration depth 
A' saturates to the system psize L. This leads to a dissi- 
pation peak at a frequencyQ 



B 
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X 



L 2 a ff L 2 



(11) 



And for very low frequencies such that 10 <C ujl , dissipa- 
tion diminishes linearly as 

X" oc ^a,. (12) 



D. Vortices with pinning 

In the presence of bulk pinning, n p ^ 0, the relevant 
time scale is tc — tff^c/^L = ')»/'»?• Now instead of 
the finite size peak at Ui, there is a vortex dissipation 
peak when tote re 1. It is worth mentioning that this 
peak frequency, 



Kp 

uj c oc — , 

Vv 



(13) 



is independent of magnetic field. For low frequencies, 
wtc <C 1, the real part of the penetration depth tends 
to A' re Ac = \/X 2 c + X L , while the imaginary part A" re 
ujtcXg /2Xc- Thus for low frequencies the dissipation 
diminishes linearly with to as 



1 



\ 2 
A c 



2ti- X C L 



u>tc- 



(14) 



Therefore, in the presence of vortices there is a new dissi- 
pation peak, in addition to the "two fluid peak" , either at 
Wi or at ujc depending on the importance of pinning. At 
low frequencies the relevant length scale for penetration 
of the ac field is either L (in the absence of pinning) or 
Ac (if pinning is important) and the dissipation is linear 
with frequency in both cases. 
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E. Surface barrier effects 

At the surface of a type II superconductor there is 
a potential barrier that prevents the entrance (and the 
exit) of vortices. This basrier has been first calculated 
by Bean and LivingstonB9 The surface barrier inhibits 
the penetration of flux at H c i, where penetration is 
thermodinamically favorable. Instead of H c i, vortices 
start to enter at the first penetration field H p > H c \. 
Above Hp, the Bean-Livinngston surface barrier can sta- 
bilize meta s tab le states with a non-equilibrium number 
of vortices.E-3cJ The interval of external magnetic fields 
where these metastable states exist in macroscopic su- 
perconductors has been obtained in Ref. p4|JI^ , f46} It has 
been found that for fields H > H p the vortex array is sep- 
arated from the boundary surface by a vortex free region 
because of the effect of-the .surface barrier. The length of 
this vortex free regiorCJEjEil is dsB — Xl cosh" 1 (H / B) . 
The ac magnetic response has been calculated by Sonin 
and Traitocl in this case. They assume that the number 
of vortices is fixed, and they allow the size of the vortex 
free region dsB to fluctuate with the rf field h ac e~ lujt . A 
new dissipation peak is found at a frequency 



USB 



tanh 2 (d SjB /A L ) 



(15) 



which for B <C H p is losb ~ V*//- ^ or ^ ow frequencies 
lo <C u>sb, the effective penetration depth goes as 



A' 



A" 



X L coth( — — ) 
sinh 2 (^) 



(16) 



Therefore, according to this approach, the dissipation 
should go as \" oc X" / L ~ lo 1 / 2 if surface barrier ef- 
fects are important. Experimentally, dissipation maxima 
attributed to surface barriers have been measured, for 
example in Ref.^7]. In platelet samples with a perpen- 
dicular magnetic field, geometrical effects enhance the 
surface barrier giving place to a "geometrical barrier" . 
Morozov et alx3 have measured the effect of geometri- 
cal barriers in the ac response of a superconductor. In 
this case there is a dissipation maximum as a function 
of magnetic field for H > H p , which is independent of 
frequency. 



III. MODEL AND DYNAMICS 

Our numerical simulations are carried out using 
the time-dependent Ginzburg-Landau equations comple- 
mented with the appropriate Maxwell equations. In the 
zero-electric potential gauge we have HE] 

9* 



dt 



(V - iA) 2 * + (1 - T)(l - I*! 2 )*] (17) 



dA 

~dt 



where ^ and A are the order parameter and vector po- 
tential respectively and T is the temperature. Equa- 
tions ([It]) and ( |l8| ) are in their dimensionless form. 
Lengths have been scaled in units of the coherence length 
£(0), times in units of to = 47rer ?l X\ /c 2 , A in units of 
H C 2 (0)£(0) and temperatures in units of T c . r\ is equal 
to the ratio of the characteristic time to for the relax- 
ation of A and the time tcL for the relaxation of 'J: 
V = t GL /t = c 2 /(4na n K 2 D), with t GL = £ 2 /D, where 
a n is the quasiparticle conductivity and D is the electron 
diffusion constant. For superconductors with magnetic 
impurities we have Di mp = c 2 /(48irK 2 a n ), and therefore 
i] = 12 in this case. 

We solve the TDGL equations following the same pro- 
cedure as in our previous work in Ref.|29|, usiag a stan- 
dard finite difference discretization scheme. e£| We dis- 
cretize space in a rectangular mesh with spacings Ax 
and Ay in each direction. The order parameter and vec- 
tor potentials are defined at the nodes of the rectangu- 
lar mesh (r = (I.J)), and the link variables U i t j — 
exp(— inh^A^ij) (/i = x, y) are introduced in order to 
maintain the gauge invariance under discretization. We 
assume that the sample has a square shape in the x, y 
direction and it is infinite in the z direction. We ap- 
ply the magnetic field parallel to the z direction, there- 
fore the problem is reduced to two dimensions neglect- 
ing all derivatives along z. The symmetry of the prob- 
lem implies for all mesh points ^4/,j = (A x ij,A y i t j,Q) 
and Bj,j = (0,0,B zIJ ), where B zI „, = (V x A) z = 
i.d x A yUJ - dyA^.j). 

The dynamical equations must be complemented with 
the appropriate boundary conditions for both the order 
parameter and the vector potential. For the order pa- 
rameter we have used the boundary condition: 



(n#)- 



(V-^)- 1 * = 



(19) 



(1 - T)Im[**(V - iA)*] - k 2 V x V x A (18) 



usually known as the superconductor-insulator boundary 
condition (S-I) because it implies that the perpendicular 
component of the superconducting current is equal to 
zero at the surface (J^ = 0). This boundary condition 
also minimizes the free energy at the sample surface. 

The ac magnetic field is introduced in the simulation 
through the boundary condition for the vector potential 
We consider the case where the ac magnetic field 
h ac is parallel to the dc-componcnt H^, and both fields 
are in the z direction: 

B z \s = (V x A) z | s = H dc + h ac cos(ujt) ; 

this expression is evaluated at the sample surface. We 
study the response of a superconductor in the linear 
regime, i.e. h ac <c Hd c - The time dependence of the 
sample magnetization can be obtained from the magnetic 
induction averaged over the sample (B(t)) through the 
relation: 

4ttM(£) = (B(t)) - H(t) 
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The ac magnetic susceptibilities are obtained from the 
Fourier transform of M(t): 

1 ( 2lx 

x'(H dc ,uj) = — — / M(t) cos(arf)d(arf) 

1 

x"(H dc ,co) = —— M{t)sm(ut)d(u;t) . (20) 

7Tft a c Jo 

In what follows, we have solved the TDGL equations 
numerically for a type II superconductor with k = 2, 
mean field temperature T = 0.5 and the parameter r\ = 
12. We used a spatial discretization of Aa; = Ay = 0.5£ 
and, in order to make efficient calculations, we have cho- 
sen adequately the time step with values At < 0.015to. 

IV. THE AC SUSCEPTIBILITY IN 
MESOSCOPIC SUPERCONDUCTORS 

In a mesoscopic sample the behavior of the ac mag- 
netic response in the presence of a dc magnetic field can 
be different than in the macroscopic case. On one hand, 
the models typically used for macroscopic samples (see 
Sec. II) , assume a high density of vortices (giving an .al- 
most uniform magnetic field profile inside the sampled), 
applied magnetic fields such that H c \ <C H <C H C 2, and 
very large (semiinfinite) samples. On the other hand, in 
mesoscopic samples there is a small number of vortices, 
surface barriers can not be neglected and the finite size 
of the sample is important. Here we perform simulations 
of squares samples of sizes from 5A x 5A to 20A x 20A 
using the TDGL equations. 

We start by showing in Fig. 1 a typical high frequency 
ac magnetic response of a small sample. In Fig. 1(a) we 
plot the driving field H(t) = H.& c + h ac cos cot and the 
corresponding magnetization M(t) of the sample. We 
show the response of the system after a long equilibra- 
tion time (typically around 100 periods of the external 
ac field). It can be observed that the magnetization has 
a sinusoidal behavior, which is a signal that we are in the 
linear regime and also that the sample is following the 
external perturbation. It can also be observed that there 
is a phase shift between the magnetization and the ex- 
ternal field clue to the presence of dissipation. Fig. 1(b) 
shows the magnetization loop obtained from the signals 
of Fig. 1 (a) . The area inside the loop is proportional to 
the time average of the energy dissipated in the sample 
and it is proportional to the imaginary part of the ac 
susceptibility x" ■ The appearance of dissipation changes 
the relationship between M ac (t) and H ac (t) and there- 
fore the real part of the ac susceptibility, x'{Hdc, oj), will 
depend on the frequency of the ac signal as well as on 
the bias field Hd c - 

A. Frequency dependence 

In order to understand the different dynamical regimes 
of the ac response, we study first the dependence of 




FIG. 1: Response of a superconductor to a high frequency 
magnetic field, (a) Time variation of the driving field and the 
sample magnetization and (b) the magnetization loop that 
shows the presence of dissipation inside the sample. (Sample 
size: 20A x 20A, H dc = 0.38# c2 (T), ut = 0.028). 



X = x' + with frequency lu for different values of 
the dc magnetic field Hd c - In Fig. 2 and Fig. 3 we show 
the imaginary and the real part of the susceptibility, re- 
spectively. These curves were obtained for a sample of 
20A x 20A. 

Fig. 2(a) shows the low frequency behavior of x" f° r 
several values of the magnetic field. Fig. 2(b) shows the 
full x"( w ) curves obtained for fields such that there are 
no vortices within the sample, i.e. in the Meissner state. 
For fields above the first penetration field H p , vortices 
enter into the sample. The full frequency dependence of 
x" (<^) in the presence of vortices, i.e. for fields Hd c > H p , 
is shown in Fig. 2(c). 

The simplest case to understand is the behavior when 
there is no external field applied, Hd c — 0. In this case 
the complex susceptibility agrees very well, both quali- 
tatively and quantitatively, with the two fluid model de- 
scribed in Sec. IIB. For low frequencies, x" nas a linear 
dependence with ui as given by Eq. ([l0]) . There is a max- 
imum in x"(w, 0) at u) = v Q = 1/fo- At high frequencies 
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FIG. 2: Frequency dependence of the imaginary component 
of the ac susceptibility x" at different fields Hdc, for samples 
of size 20A x 20A. (a) Low frequency behavior for fields in the 
Meissner state (black symbols, Hdc/ H C 2(T) = 0,0.30) and in 
the mixed state (open symbols, H dc /H c2 (T) = 0.40, 1.0). The 
continuous lines correspond to a linear dependence \" ~ u). 
(b) Frequency dependence in the Meissner state (number of 
vortices: N v — 0) The dashed line corresponds to a depen- 
dence x" ~ w _1 . (c) Frequency dependence in the mixed 
state (N v ^ 0). The frequency is normalized by v a = 1/to- 

u> ^> vq the dissipation is due to the normal electrons and 
follows the expected dependence x" ~ u~ x / 2 as given by 
Eq. (||), corresponding to the normal state skin depth 
effect. 

For small magnetic fields within the Meissner state, 
< Hdc "C H p , the behavior is qualitatively similar to 
the Hdc — case, with a peak at u = v§. The main 
difference is that the slope of the low uo linear depen- 
dence increases with magnetic field, as can be seen in 
Fig. 2 (a). According to Eq. (jl(]), the linear slope for low 
frequencies is proportional to Al- This implies that the 
London penetration depth should be dependent on mag- 
netic field: Xl(H). Also, Fig. 2(b) shows that the magni- 
tude of x" (w) at any given tu increases for increasing Hd c - 
These two effects are easily understood since the Meiss- 
ner screening currents deplete the magnitude of the order 



FIG. 3: (a) Frequency dependence of the real component of 
the ac susceptibility x a t different values of Hdc (black sym- 
bols: Meissner state, open symbols: mixed state) and (b) 
effective ac penetration depth \ ac as a function of frequency. 



parameter at the surface, which leads to a larger static 
penetration depth Xl{H) (see also Sec.IVB) and to an 
increase in the normal electron dissipation. The most in- 
teresting result is that for Hd c near but still below H p a 
second dissipative maximum appears in %"(a;). This sec- 
ond peak is at a frequency lu p two orders of magnitude 
below v a . 

For Hdc > Hp the vortices enter into the sample. In 
this case we still find a low frequency linear dependence 
X" ~ oj, as shown in Fig. 2(a). Only for large fields 
close to H C 2 a small departure from linear dependence 
is observed, x" ~ a - ,Q i with a « 0.8. The behavior for 
very large frequencies w > is the same for all mag- 
netic fields, since it corresponds to the normal state skin 
depth effect with x" ~ ^ _1 ^ 2 , as mentioned before. As 
can be observed in Fig. 2(c), the main characteristic of 
the mixed state is that there are two dissipation maxima 
as a function of frequency in x"(o;). The high frequency 
maximum corresponds to the "two fluid peak" at uj w vq. 
The low frequency maximum that appears at a frequency 
cl> p <^ vo, is a continuation of the maximum observed 
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when Hdc < H p . The peak frequency uj p shows a weak 
dependence with magnetic field, uj p first decreaseas and 
then increases with increasing Hd c - For fields Hd c > H c3 
the sample is completely in the normal state and there is 
a single maximum in the x" vs 10 curve. It corresponds 
to the frequency at which the normal state skin deph 8 n 
equals the system size L. 

Let us now analize the frequency dependence of the 
real part of the susceptiibility, x'(w), which is shown in 
Fig. 3(a) for some values of Hdc- From x'(u>) it is pos- 
sible to extract the real part of the effective ac penetra- 
tion depth X ac = X'(oj) using Eq. (||). The length \ ac 
represents the length scale for the penetracion of the ac 
component of the magnetic field and its frequency de- 
pendence is shown in Fig. 3(b). In the absence of a bias 
field {H dc = 0), x'( w ) is we U described by the two fluid 
model, as given in Eqs. (|9|) , (|l(]) . The sample changes 
from a London screening for u> <C v a , with \ ac w A^, to 
a perfect diamagnetic state for u v a , with A ac «s 0. At 
low frequencies x' is larger than — l/4n due to the pene- 
tration of the field in a region of size \l near the sample 
surface, x' ~ ( — l/47r)(l— 4Xl/L). At high frequency the 
ac penetration depth is A ac < Ai and x'(a;) ~ —1/4%. 
The opposite case is for Hd c = 1.9H C 2(T) > H cS with 
the sample mostly in the normal state (at this field the 
surface superconductivity disappears and thexe is super- 
conductivity only at the corner of the sampleuj) . As can 
be expected, in this case the system behaves as a nor- 
mal metal and is dominated by the skin-depth effect with 
A Q c ~ S n /2. For low frequencies, the skin-depth 5 n equals 
the system size and thus \' « 0. For magnetic fields 
within the Meissner state, < Hdc *C H p , we find that 
for low frequencies, wCvoj the ac penetration depth sat- 
urates to the field dependent London penetration depth 
\l(H), which increseas with H. The most interesting 
case is for Hd c < H p , which is shown in Fig. 3. In this 
situation, one can distinguish three frequency regimes, as 
it was found before in the dissipative part, At the 

frequencies for which there are maxima in the dissipation, 
Vq and u) p , there is a rapid change in the value of x'( w ) 
(or in X ac (uj)). For lo 3> Vq the normal fluid behavior 
of Eq. (H) is followed, as expected for all values of Hd c . 
For the frequency range io v < lo < i/q, there is a shoulder 
in x'i 1 - ) (and A ac (o>)), which basically corresponds to an 
ac penetration depth of the order of the expected value 
for Xl(H) at that field. At low frequencies, lu < u p , the 
ac penetration depth saturates to a length i p which is 
smaller than the system size and larger than A/, (if). 

In summary, the most interesting feature observed in 
x(to) is the appearance of a dissipation maximum at fre- 
quencies lu p (H) <C vq in the presence of a magnetic field. 
A dissipation maximum in macroscopic samples is ex- 
pected when H ^ 0, as we discussed in Sec. II. It can 
be caused by several reasons: (i) finite size, (ii) surface 
barriers or (iii) bulk pinning. However, the behavior pre- 
dicted in any of these cases is not followed in the meso- 
scopic samples studied. Let us discuss each of the possi- 
bilities, (i) Finite size: in a small sample this is the first 
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FIG. 4: Dependence with the TDGL time scale tch = r/to. 
x" vs. lo curves for small mesoscopic samples 5A x 5A in the 
Meissner state (N v = 0) near the penetration field, H < H p , 
Hdc = 0.35.ff C 2(T). Inset shows the dependence with 77 for 
H dc = 0. 



effect to analyze. The characteristic frequency for finite 
size effects, ljl, should increase linearly with B and de- 
crease with system size, as given by Eq. (0). However, 
ujp has a weak and non-monotonous dependence with 
magnetic field, as seen in Fig. 2(c). In Fig. 4 we show 
x"(u;) in a smaller sample with L — 5A (and i] = 12) 
for Hdc ^ H p . We see that the maximum appears at 
the same frequency uj p as in the L = 20A case for the 
same magnetic field [Fig. 2(b)]. Therefore lo p has no size 
dependence. Furthermore, in Fig. 3 we see that when 
lo — > the sample is still diamagnetic and A ac < L (a 
finite size effect would have x'( w ^ 0) w and A ac > L). 
(ii) Surface barrier: in mesoscopic samples, surface bar- 
riers are very important.Ej In macroscopic samples, the 
effect of surface barriers, assuming a large number of vor- 
tices, gives a low frequency dependence for dissipation as 
X" ~ uj a with a = 1/2, see Eq. @ and Ref. |. How- 
ever, we find a linear frequency dependence in most cases 
as shown in the inset of Fig. 2(a), and even when there 
is a departure from linearity, it is with an exponent of 
a w 0.8 > 0.5. Also the characteristic frequency ojsb 
of Eq. ( |i~5| ) has a strong dependence on magnetic field, 
which is not observed here, (iii) Pinning: since there is 
no disorder in the model, no bulk pinning is expected. 
However, in a very small sample one may argue that, 
given the smallness of the system, the sample itself acts 
as a confinement potential for vortices. In this case, k p 
in Eq. (0) would represent the vortex elastic response 
for small oscillations within this confinement potential. 
The characteristic frequency for pinning effects, u>c, is 
independent of magnetic field in macroscopic samples, 
according to Eq. (|l3|). The weak field dependence of lu p 
may be still consistent with this result (the effective n p 
may depend on the number and distribution of vortices 
inside the sample in this hypothetic case of "confinement 
potential"). Also the fact that A ac < L for w « w p in 
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FIG. 5: History dependence of x vs - w curves for Hdc near 
the penetration field H p (Hdc = 0.38H C 2(T)). Black triangles: 
decreasing frequency from a high ui state without vortices 
(Nv = 0); the branch for u) above the jump has N v = 
and below the jump N v — 12. Black symbols: increasing 
frequency from a low ui state with N v — 12. Sample size is 
20A x 20A. The inset shows the absence of history dependence 
for a,H d c = 0-GHc2(T). 



tices into the sample, leads to hystory dependent effects. 
The curves in Fig. 2 were generated increasing the field 
Hdc at low frequency and then measuring x" (<^0 for in- 
creasing u>. In this way, the curves reported in Fig. 2(c) 
have the same number of vortices at all frequencies for 
a given Hdc- If we follow the opposite procedure, i.e., 
we increase Hd c at high frequency and then measure 
X"(w) for decreasing oj, the curves can show discontinu- 
ous jumps at certain frequencies values due to vortex pen- 
etration. This occurs because increasing the field at high 
frequency can generate a metaestable state that can be- 
come unstable at low frequencies. In Fig. 5 we show this 
case, for a field near the first penetration field Hd c ~ -Hp- 
Depending on the history, the high frequency branch can 
have either N v = 12 vortices or N v = 0. In the case when 
the frequency is decreased from a high frequency state, 
there is a large jump in at a frequency near u> p , 

from a metastable branch with N v = for u> > lo p to a 
branch with N v = 12 for uj < uj p . In contrast, we see in 
the inset of Fig. 5 that for a field Hdc ^ H p the curves of 
x"(w) are independent of history. 



B. Magnetic field dependence 



Fig. 3(b) is consistent with a finite "Campbell's penetra- 
tion depth" Ac (and then t p = Ac)- However, we also 
have to discard this possibility, since we find that there 
is a dissipation maximum at lo p even when there are no 
vortices in the sample. As shown in Fig. 2(b) the dissipa- 
tion peak already appears for magnetic fields just below 
the penetration field H p . To be more precise, we have 
explicitly calculated the topological number of total vor- 
ticity as N v = J- §(A+^)dl. We obtain N v = at all 
times for the magnetic fields shown in Fig. 2(b). There- 
fore, all arguments based on the oscillation-. of vortices, 
as for example the Coffey and Clem modelta of Eq. (j^), 
which work very well for macroscopic samples, are not 
enough to explain the dissipation maxima observed in 
these mesoscopic samples. 

In Fig. 4 we show that oj p depends directly on the 
time scale for the relaxation of the amplitude of the or- 
der parameter, icL = V^o- We consider a magnetic field 
Hdc H p and we change the value of 77 in Eqs.(|l7|) and 
(|l8|). The "two fluid peak" is always at the same fre- 
quency vq — l/to, since to is fixed. We find that when 
decreasing 77, the frequency uj p shifts to higher values and 
increases monotonically with 1/tcL, until the rather un- 
physical case of toL < to (v < 1)> when only the two 
fluid peak at vo is observed. The dependence of u> p with 
tQL shows that variations in the amplitude of the order 
parameter are the main mechamism of this dissipation 
maximum, as we will see more clearly later in Sec. V. In 
the inset of Fig. 4 we verify that for Hd c = the behav- 
ior is independent of tcL, since this case can be simply 
described with the two fluid model. 

The appearance of the dissipation maximum for a mag- 
netic field below H p , just before the penetration of vor- 



Experimentally, the ac frequency is fixed and the dc 
magnetic field or the temperature can be varied. In 
this section we will study the ac magnetic response of 
mesoscopic superconductors varying the bias dc magnetic 
field. We will show results at different frequencies and for 
different sample sizes. 

In the absence of the ac field, the dc magnetic behavior 
of mesoscopic superconductops,|is, different from the con- 
tinuos macroscopic behavior.! 16 ^ 9 ! In mesoscopic samples 
each vortex entrance event produces discontinues jumps 
in the magnetization curve at certain Hd c values. Each 
discontinuous jump in M(H) corresponds to a sudden 
increase in the number of vortices. These jumps occur at 
succesive magnetic fields H p i = H p \, H p %, H P 3, . . .O In 
the regions of H p i < H < H p ^ i+ ^ the number of vortices 
N v is constant. The only penetration events occur at 
H p i, when vortices enters into the sample. In the region 
of constant vorticity, H pl < H < H p u+i\ , one may think 
that no vortices enter the sample because of the effect of 
the surface barrier. 

Fig. 6 (a) shows the behavior of x" (Hdc) obtained at a 
fixed frequency ui — 0.004i/ o < uj p in a mesoscopic sample 
of size 10A x 10A. At small Hd c , in the Meissner state, 
X" increases continuously with increasing Hd c - The pres- 
ence of the dc magnetic field in the Meissner state induces 
static supercurrents which deplete the order parameter at 
the boundary. As explained before in Sec. IVA, the de- 
pletion of the order parameter leads to an increase of dis- 
sipation for increasing Hd c - When the dc magnetic field 
is increased above the first penetration field H p = H p \, 
the first vortices enter into the sample. Fig. 6(a) shows 
that there is a discontinuous jump in x" with a decrease 
of dissipation just at H p . Two states are possible exactly 
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FIG. 6: Hdc dependence of x" obtained at a fixed frequency. 
Small mesoscopic samples of size 10A x 10A showing discon- 
tinues jumps: (a) u)to = 0.004 and (b) uto = 0.06. (c) Large 
samples of size 40A x 40A showing the usual continuos behav- 
ior, Loto — 0.06. 



at H p : one state without vortices and high dissipation, 
and one sate with vortices inside the sample and low dissi- 
pation. This is in agreement with the history dependence 
observed in Fig. 5. Further jumps in the dissipation curve 
X" (Hdc) are present at the other magnetic fields for vor- 
tex penetration, H P 2, H P 3, .... The jumps are followed by 
a later continuos increase of x" (Hdc) with increasing Hd c 
while the number of vortices N v remains fixed. Fig. 6(b) 
was obtained in a sample of the same size of Fig. 6(a) but 
at a higher frequency oj = 0.06i' o > ui p . The frequencies 
used in Figs. 6(a) and 6(b) are at both sides of the dis- 
sipation maximum that appears at u p . Both curves are 
similar in their qualitative features, but show a few dif- 
ferences. At high magnetic fields there are less jumps in 
x"(Hdc) in Fig. 6(b) than in Fig. 6(a). This is because at 
high frequencies it is possible to remain in a metastable 
state with a fixed number of vortices in a wider range of 
magnetic field. At the same time, the higher frequency 
of Fig. 6(b) produces a decrease in the amplitude of the 
jumps. The behavior of x"(Hd c ) of Figs. 6(a) and 6(b) 



FIG. 7: Hdc dependence of x obtained at a fixed frequency: 
(a) and (b) for small samples and (c) for large samples. The 
parameters are the same as used in Fig. 6. 



is different from the behavior of bulk samples. For ex- 
ample, Fig. 6(c) shows x" (Hdc) at w = 0.06^0 obtained 
in a "large" sample of 40A x 40A. In this case there are 
no discontinuous jumps and the curve is almost contin- 
uous. In macroscopic samples the entrance of vortices 
for Hdc > H p increases the ac losses in a continuos way. 
In the macroscopic case vortices play a fundamental role 
in the ac losses, and therefore an increase on the num- 
ber of vortices increases the dissipation. In contrast, the 
results of Figs. 6(a) and 6(b) show the opposite behav- 
ior in mesoscopic samples: there is a sudden decrease of 
the ac losses at the magnetic fields that correspond to an 
increase in the number of vortices. 

Fig. 7 shows the real part of the ac susceptibility for 
increasing Hdc- We use the same parameters as in Fig. 6. 
The curves of x'(Hdc) also have a typical mesoscopic be- 
havior with discontinuities at each vortex penetration 
field. The increase of frequency from ui = 0.004i^ o in 
Fig. 7(a) to uj — 0.06^0 in Fig. 7(b) shows similar dif- 
ferences as observed in x"(Hdc)- It can also be noted in 
Fig. 7(a) that the differential susceptibility has a para- 
magnetic behavior (%' > 0) near each vortex penetra- 
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FIG. 8: Change in the ac penetration depth (AA ac ) in a meso- 
scopic sample (squares, 10A x 10A) and a macroscopic sample 
(circles, 40A x 40A) for uito = 0.06. The inset shows a fit to a 
quadratic dependence AA ac ~ H 2 (continuos lines) for fields 
below Hp. 

tion field. This can be expected from the behavipc of 
the dc magnetization curves of mesoscopic samples Jl3 In 
mesoscopic samples, for H dc below each vortex penetra- 
tion field a magnetization maximum appears and there- 
fore a region near H p i where the magnetization increases 
at increasing field. This shows in the low frequency ac 
response of mesoscopic samples as a differential param- 
agnetic behavior. Fig. 7(c) shows a continuos macro- 
scopic dependence of x'(H dc ) in a larger sample of size 
40A x 40A. At the frequencies used in Fig. 7(c) we can 
see a diamagnetic behavior even in the normal state for 
fields greater than 2H C 2(T). This is simply because at 
high frequencies there is magnetic screening in the nor- 
mal state due to the skin depth effect. 

From Fig. 7 it is possible to obtain the H dc dependence 
of the London penetration depth. This dependence has 
been a subject of interest recently because of its relation 
with the simmetry of the order parameter. In conven- 
tional s-wave superconductors the following dependence 
is expected: 

AX(H dc , T)/A(0, T) = f3(T)(H dc /H c2 (T)) 2 (21) 

where AX(H dC7 T) = X{H dc ,T) - A(0,T). It is assumed 
that the sample is in the Meissner state and in the lin- 
ear regime, i.e. \H dc \ > \h ac \. Eq. |jj|) was origi- 
nally obtained by Ginzburg and LandauEj and also by 
BajdeenEJ an it has been observed experimentally in type 
iBB and conventional type II superconductors!!! The be- 
havior of unconventional superconductors is expected to 
be differentia. r-LoJlie case of a d-wave superconductor 
Yip and Saulsti(E3 proposed a linear field dependence 
AX(H dc ) oc H dc . Here we use the phenomenological 
TDGL equation in its s-wave form. Therefore, in the 
Meissner state we expect a result similar to Eq.(^Tj) for 
macroscopic samples. The results are shown in Fig. 8. 
Both the macroscopic and the mesoscopic sample show 



a quadratic dependence of AX(H dc ) as a function of H dc 
in the Meissner state. The vortex penetration produces 
a sudden decrease of AA in mesoscopic samples, whereas 
at a fixed number of vortices AA increases continuously 
for increasing H dc . For macroscopic samples the entrance 
of vortices at the first penetration field H p does not pro- 
duce discontinuity in AX(H dc ) but a change of slope as 
can be observed in the plot. For H dc > H p the depen- 
dence changes from quadratic to approximately linear in 
H dc . 



V. NASCENT VORTICES EFFECT 

From the results of the previous section, it is clear now 
that the ac magnetic response of a mesoscopic supercon- 
ductor is different from the macroscopic behavior and it 
can not be explained with the known models. An inter- 
esting result is the significant decrease of dissipation each 
time the number of vortices increases. To understand 
this effect we have analyzed in detail the time variation 
of the order parameter and the magnetization within the 
sample. 

Fig. 9 shows the time variation of the order parame- 
ter at different applied magnetic fields in a small sample. 
Each curve represents the amplitude of the order param- 
eter along the x direction, perpendicular to a face of the 
square sample, obtained at a given time. To simplify 
the analysis the cut was done just at the center of the 
face at y — L/2 (which is the preferred point for vortex 
penetration in a square sample). The order parameter is 
obtained at different times within a single period of the ac 
magnetic field. In particular, we show the behavior of the 
order parameter just before and after the two first values 
of the penetration fields, H v \ and H P 2. At H dc = the 
order parameter 1^1 is equal to unity in the entire sample 
and the dissipation is small. In the Meissner state but for 
H dc < H p i, Fig. 9 (a) shows that the magnitude of the 
order parameter near the boundary, and at the center of 
the face, has large oscillations around a small value. Fig. 
9 (b) shows the time evolution of order parameter just 
after the first vortex entrance event, H dc > H p %. We see 
that the entrance of vortices increases the order param- 
eter at the boundary, where now oscillates around a 
value much higher than before. One can also see that the 
positions of the vortices are fixed inside the sample and 
do not oscilllate in the presence of the ac signal. Recently, 
in an extension of the idea of the Bean Livingston barrier 
to mesoscopic superconductos, we have shown that the 
vortex entrance produces a reinforcement of the surface 
barrier, due to the interaction between the vortices that 
are trying to entex_.the sample and the vortices already 
inside the samplers The presence of this interaction de- 
creases the surface currents. Thus, the main effect after 
the entrance of vortices inside a mesoscopic supercon- 
ductor is that the Meissner currents decrease and the or- 
der parameter increases at the boundary. Therefore the 
normal state electrons decrease and then the dissipation 
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FIG. 9: Time variation of the order parameter in a sam- 
ple of size 10A x 10A at different applied magnetic fields and 
for u)to = 0.06. The figures show transversals cuts taken at 
the center of one of the faces for different times, (a) Hdc = 
0.4216# o2 (T) < H p i, (b) H ic = 0.4217# c2 (T) > H p i, (c) 
H dc = 0.568H c2 (T) < Hp2, (d) H dc = 0.569# c2 (T) > H p2 . 

decreases. 

Figs. 9(c) and 9(d) show similar behavior for Hd c < 
H p 2 and Hd c ^ H v i- After increasing Hd c at a constant 
number of vortices, the order parameter at the boundary 
has decreased, as seen in Fig. 9(c). The decrease of the 
order parameter increases the dissipation as was observed 
before. If new vortices enter again into the sample, the 
order parameter increases at the boundary, see Fig. 9(d) 
and therefore the dissipation decreases again. This pro- 
cess is repeated in the other penetration fields H pil . On 
the contrary, the entrance of vortices in a large sample do 
not produce an appreciable change in the amplitude of 
the order parameter at the boundary before and after the 
penetration event, and therefore the dissipation shows a 
continuous behavior [as seen in Fig. 6(c)]. 

A careful analysis of the entire profile of the order pa- 
rameter for magnetic fields below but near to H p sug- 
gest that the depletion of the order parameter near the 



boundary may be attributed to "vortices" that are about 
to enter into the sample. These "nascent vortices" are lo- 
cated outside of the sample but near the surface. Even 
more, for higher values oih ac (beyond linear response) we 
find that these "nascent vortices" oscillate near the sur- 
face (i.e. that the profile of the depletion of |\E'| near the 
surface oscillates with H(t)). According to this scenario, 
the length scale t v observed in Fig. 3(b) for low frequen- 
cies may correspond to a characteristic length scale for 
the nucleation of vortices from the boundary. (For large 
fields B ^> H p , the length l v could be related to the 
size of the vortex free region dsB of macroscopic sam- 
ples). Years ago, Walton and RosemblunS suggested 
the idea of nascent vortices as a possible source of the 
high-frequency losses in superconductors. This idea had, 
further support from the analytical work of Kramero 
who found static solutions of the Ginzburg-Landau equa- 
tions which corresponded to vortices nucleated near the 
sample. The nascent vortices can oscillate under an ac 
field. This source of dissipation is not generally taken 
into account in macroscopic samples because in this case 
is small compared with the dissipation due to the motion 
of vortices inside the sample. However in a mesoscopic 
superconductor, with a few vortices confined inside the 
sample, this source of dissipation can be a significative 
part of the dissipation observed, as we have found here. 



VI. SUMMARY AND CONCLUSION 

We have studied the response of mesoscopic supercon- 
ductors to an ac magnetic field. In mesoscopic super- 
conductors vortices are confined inside the sample by the 
surface barrier and their dynamics play a secondary role 
in the dissipation. When increasing Hd c , discontinuities 
in x' an d x" appear at each vortex penetration event. 
The dissipation is maximum before the penetration of 
vortices due to the effect of "nascent vortices" , which lead 
to large oscillations of the amplitude of the order param- 
eter at the boundary. After the penetration of vortices, 
the normal currents decreases and therefore the dissipa- 
tion decreases, leading to a sudden suppression of the ac 
losses induced by the vortex entrance. 
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